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The probability-generating function of the number of electron-positron pairs produced in a uniform 
electric field is constructed. The mean and variance of the numbers of pairs are calculated, and 
analytical expressions for the probability of low numbers of electron-positron pairs are given. A 
recursive formula is derived for evaluating the probability of any number of pairs. In electric fields 
of supercritical strength \eE\ > 7rm 2 /ln2, where e is the electron charge, E is the electric field, and 
m is the electron mass, a branch-point singularity of the probability-generating function penetrates 
the unit circle \z\ = 1, which leads to the asymptotic divergence of the cumulative probability. 
This divergence indicates a failure of the continuum limit approximation. In the continuum limit 
and for any field strength, the positive definiteness of the probability is violated in the tail of the 
distribution. Analyticity, convergence, and positive definiteness are restored upon the summation 
over discrete levels of electrons in the normalization volume. Numerical examples illustrating the 
field strength dependence of the asymptotic behavior of the probability distribution are presented. 



PACS numbers: 11.15.Kc, 12.20.Ds 

I. INTRODUCTION 

Klein's paradox [l| is an excellent illustration of the 
problems associated with the existence of solutions with 
negative energy in the Dirac equation. Since the formu- 
lation of the paradox, it has attracted considerable at- 
tention. The physical meaning of the apparently anoma- 
lous behavior of the solutions of the Dirac equation has 
gradually been clarified, and the complete solution to the 
paradox has finally been obtained 0-@ • The effect is re- 
lated to the production of electron-positron pairs in an 
external electric field or, equivalently, with the possibility 
of the tunneling of electrons from the region containing 
solutions with negative energy to the region containing 
solutions with positive energy. 

Verification of the Schwinger mechanism of electron- 
positron pair production |5l4lOl| is still an open experi- 
mental challenge, due to the smallness of the effect for the 
electric fields that are available in a modern laboratory. 
In recent years, the mechanism of pair production based 
on multiple soft-photon absorption has been discussed. 
The combination of a strong, slowly varying field and a 
weak but dynamical field lowers the threshold for pair 
production, which provides a possible method to create 
electron-positron pairs in experiments with ultra-high- 
intensity lasers [11] . Various other possibilities to create 
electron-positron pairs in laser beams are also discussed 
[12l 4l5| . Latest developments in this field are reviewed in 
Ref. 

The probabilities of creating one, two, or more 
electron-positron pairs that contain detailed information 
about the pair-production mechanism fundamentally al- 
low for a direct experimental measurement. In this paper, 
we construct the probability distribution of the electron- 



positron pairs produced in a uniform static electric field. 
In Sec. II, the probability distribution is calculated on 
the basis of the probability-generating function (PGF), 
and in Sec. Ill, the probability distribution is calculated 
with the help of recursion relations. As shown in Sec. IV, 
the continuum limit approximation is restricted by weak 
electric fields. In strong electric fields, the analytic prop- 
erties of the PGF are violated and the cumulative prob- 
ability is divergent. We show that the summation over 
discrete modes restores both the analyticity of the PGF 
and the convergence of the series. The continuum limit 
approximation also leads to a breakdown of the positive 
definiteness of the probability in the tail of the distri- 
bution. This shortcoming is also overcome in a discrete 
scheme. The main features of the asymptotic behavior of 
the probability distribution are illustrated with numeri- 
cal examples. Final conclusions are collected in Sec. V. 



II. THE PROBABILITY-GENERATING 
FUNCTION 

In the standard formulation of the problem, the electric 
field is switched on adiabatically for a time T, and then 
the field is adiabatically switched off. The probability 
that a vacuum remains a vacuum Q and the average 
number of produced pairs @ are of interest. We are also 
interested in the probability distribution of the number 
of produced pairs. 

Electrons with negative energy moving against the 
electric field lines are slowed down and stopped, and 
they either tunnel with a probability of w a into positive- 
energy states or turn back with a probability of 1 — w a , 
in which case they are further accelerated along the field 
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lines. We denote by a = (p, a) the quantum numbers of 
electrons, where p is the initial momentum, and a is the 
electron spin projection. In a uniform electric field, the 
semiclassical expression for the probability of tunneling 
appears to be accurate: 



exp 



7r(m 2 + p 2 ± ) 



(1) 



where e is the electron charge, E is the electric field 
strength, and p^ is the electron momentum orthogonal 
to the field. 

Thus, we have an example of Bernoulli trials: the tun- 
neling of electrons occurs with a probability of w a for 
each electron, and the reflection of electrons occurs with 
a probably of 1 — w a for each electron. 

The vacuum-to-vacuum transition probability equals 



Po = - Wa). 



In electric fields, any number of electron-positron pairs 
can be produced. The sum of all of the probabilities is 
equal to unity: 



U(i - w a ) + n (i - 

a a 

+ 27 J2 WaW P II (1 - «>t) + • • • = L ( 2 ) 

Whenever the series converges, this algebraic identity is 
independent of w a . The proof of this fact can easily be 
made using a representation of the series in the form of 
a product, namely, 



PO JJ(l + ^a), 



(3) 



where v a — w a /(l — w a ). The sequential terms of the 
left-hand side of Eq. ([2]) have the meaning of probability 
that there are n = 0, 1, 2, ... electron-positron pairs. 
Consider a more general expression: 



<p{z) 



Po 



(4) 



In probability theory and statistics, the function (p(z) is 
known as the PGF. A stochastic process to which it cor- 
responds is characterized by the probability distribution 



I d n 

P " = d^ (z)l ^°- 



(5) 



In our case, each of the quantities p n describes the prob- 
ability of there being n electron-positron pairs. 

A PGF has many useful properties. In particular, 
ip(z) = Fi[z n ] is the expectation value of z n ; the equa- 
tion ip(l) — 1 is equivalent to Eq. ([5); <p(Q) = p is 
the vacuum-to-vacuum transition probability; the aver- 
age number of electron-positron pairs is E[n] = ¥>'(!); the 



variance equals V[n] = E[n 2 ] — E[n] 2 = ip"(l) + y>'(l) — 
^'(l) 2 . 

The logarithm of <p(z) includes a summation over the 
indices a. During the period of observation T, the rel- 
ativistic electrons with longitudinal momenta < pu < 
\eE\T slow down and then stop their movement along 
the field lines. Once they have stopped, these electrons 
acquire a chance to tunnel. In the continuum limit ap- 
proximation, computing the discrete sum is replaced by 
integrating over the phase space: 



2Vd 3 p 



VT\eE\ 



2d 2 p A 



(6) 



where V is the normalization volume. The probability of 
tunneling does not depend on pu , so this component was 
integrated out. 

Using well-known arguments 043' we obtain 



tp(z) = exp (-7Li 2 ((l - z)0) , 
where Li2(t) is the dilogarithm 



(7) 



U 2 (t) = -J 



ln(l - 0) 



de 



and 



7 = ^^, ^ = cx P (-— -). 
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The value of 7 is the effective number of negative-energy 
electrons that reach a turning point during time T and 
then obtain a chance to tunnel. To replace the discrete 
sum with an integral, one has to require 7 > 1. As we 
shall see later, this condition, although necessary, is not 
sufficient. HU 

The vacuum-to-vacuum transition probability and the 
probabilities of producing one, two, or three electron- 
positron pairs can easily be found by means of Eq. ([5} : 



Po 
Pi 

P2 = 
P3 = 



exp(-7Li 2 (£)) , 
7^Po, 



7 



7 i 2 



L - 



Po, 



(l-o 2 



-2 + 2<e + 3 7 <e 



L 



Po, 



(8) 
(9) 

(10) 



(11) 



where L = — ln(l — £). The zero-order term po is in 
agreement with the result of Schwinger Q . The average 
number of pairs, which was first found by Nikishov [6j], 
equals 



E[n] = 7 £, 



(12) 
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FIG. 1: (color online) The complex z-plane of the PGF for 
electric fields of subcritical strength. The branch-point singu- 
larity of the PGF at z = —r^ is outside the unit circle \z\ = 1. 
The contour C shows the path of integration in Eq. (| ltjp : the 
integrand also has a pole singularity at z — 0. 

while the variance in the pair number is given by 

V[n] = i 7 £(2 - £). (13) 

The leading terms of p n /po ~ ( , yL) n /n\ induce the 
Poisson distribution. They describe the n uncorrelated 
pairs and, in violation of fermion statistics, they count 
configurations with two or more negative-energy elec- 
trons in the same state. The lower-order terms in 7 re- 
move these unphysical contributions. For example, the 
probability of there being two pairs, P2, is given by 

^ WpPo = 2^ "<*) 2 Po ~ 5 ^2( v ") 2 Po- ( 14 ) 

a a 

The first positive term equals po(^L) 2 /2l; the second neg- 
ative term has the order 0(7) (when it is divided by 
po). The left-hand side is positive, so the second term is 
smaller than the first one. To avoid violating this con- 
dition, the value of £ should not be especially close to 
unity: 

l-£>l/7- (15) 

The expansion in the powers of 7 in the expressions for 
Pn/po arises due to fermion statistics, and this expansion 
has a combinatorial meaning. The production of uncor- 
related pairs is described by the Poisson distribution with 
the mean value E[n] = jL. The correlation effects show 
up in the deviation of Li2(£) from L in the exponent of po, 
in the deviation of £ from L in the mean value, etc. The 
pre-exponential factor of p\ has a purely Poisson form, 
since in the one-pair production the correlations are not 
important. 



FIG. 2: (color online) The complex z-plane of the PGF for 
electric fields of supercritical strength. The branch-point sin- 
gularity of the PGF at z = —r^ penetrates the unit circle 
\z\ = 1, which destroys the convergence of the cumulative 
probability. The contour C shows the path of integration in 
Eq. CSJ. 

III. RECURSION RELATIONS 

The calculation of p n for large values of n can be per- 
formed recursively using the method employed in Ref. 




The function Li2(i) is analytic in the complex i-plane, 
except for the logarithmic branch point at t = 1. The 
function (p(z), which is expressed in terms of Li2(i), also 
has a logarithmic branch point at z = —r^ = 1 — l/£. 
For £ < 1, <p(z) is analytic in the neighborhood of z = 0; 
therefore, we can represent p n by the following contour 
integral: 

1 f dz 

Pn = ^i! ^ {z) - (16) 

c 

The contour C encompasses the point z = 0, as shown in 

Figs. Q] and [2j As long as ip(z) is regular at z = 0, p n = 

for n = —1, —2, ..... 

To obtain the recursion, we carry out integration by 

parts in Eq. ([r6]) : 

1 f 1 

p n = 7T~- / — r<£ (*)dz. 
2m J nz n 

c 

This expression also follows directly from the definition 
of PGF. The first derivative of ip(z) has a simple form: 

ip'{z) = In (1- ¥>(*). (17) 

1 — z 

By decomposing the logarithm in front of f(z) first in 
a power series of 1 — z, then in a power series of z, we 
obtain 

fe=l 1=0 v ' 
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One can carry out the summation explicitly: 

k=l+l x ' 

= (-i) l f^F(i + hi + i;i + 2;0, (is) 

where F(a, 6; c; z) is the hypergeometric function. For 
/ = 0, this expression gives go = L; for / > 1 wc obtain 

m=l 4 

Finally, one finds 

n-l 

Pn = - VWn-l-1- (19) 

n 

This equation gives the recursive definition of the proba- 
bility distribution p n of the random variable, n, which 
represents the number of electron-positron pairs pro- 
duced in a uniform electric field. 
Because 

oo 

^2 91 =t, 

l=Q 

we reproduce Eq. (|12j) after the summation of np n over 
the pair number. 

IV. ASYMPTOTIC BEHAVIOR OF THE 
PROBABILITY DISTRIBUTION 

The radius of convergence of the Taylor series expan- 
sion of <f{z) at z = is determined by the distance to 
the nearest singularity. The series over z allows one to 
evaluate <p(z) at z = 1, provided = l/£ — 1 is greater 
than one. This condition is satisfied for £ < 1/2. Under 
this condition, one can rigorously prove that the sum of 
the probabilities §5$ equals one. 

The appearance of the branch-point singularity inside 
the unit circle destroys the precise probabilistic meaning 
of the PGF. The only explanation we can offer for this 
phenomenon is that the infinite product (J3J diverges in 
the continuum limit. If this divergence did not occur, 
the radius of convergence would be at least 1. Although 
product ((SJ diverges, product (j4|), which differs from ([3]) 
only by a factor of z for v a , converges for z inside the 
circle \z\ < < 1. If 1 < the first product converges; 
the second product converges inside the circle \z\ < rj. 
Thus, although the p n of equation ([5]) are analytic func- 
tions of the parameter £, and although equation ([5]) per- 
mits an analytical continuation to the region £ > 1/2, the 
physical meaning of the p n in strong fields is uncertain. 
Evidently, the applicability of equation ([5]) is limited by 
the condition £ < 1/2. 



The effective number of electrons with negative energy 
that have a chance to tunnel was previously estimated as 
~ 7. Therefore, it is expected that the p n are vanishingly 
small for n > 7. We also have a condition that the term 
( , yL) n /ri\ in the expansion oip n /po over 7 is the leading 
one, and therefore it obviously exceeds the lowest-order 
term ~ 7. The estimate of the latter gives (cf. Eqs. (j9]) 
- (TED and (HI) 

AWpq ~ (- 1 )"" 1 ~ (-l)^Vr-r 1 - ( 20 ) 

a 

The condition \Ap n \/p < (jL) n /nl is equivalent to 

n < n c = jLr^. (21) 

This inequality is fulfilled for all n < 7, so, again, we find 
> 1, which is in agreement with the previous analysis. 
For n = 2, we return to the less stringent constraint 
(|15|) . The condition > 1 also provides a term-by-term 
convergence of the sum of p n . For large n, the inequality 
pTjl is violated in any case. This violation does not lead 
to contradictions as long as the sum of the p n is saturated 
for n < 7. 

Apparently, the reason that the p n do not sum up to 
unity for < 1 is that the replacement ([6]) is not jus- 
tified. The subtraction of terms, which restores fermion 
statistics, gives rise to contributions of the type (|2"0)) . For 
n ~ 7, the crucial point is whether the magnitude of v a 
is greater than unity or less than unity. If v a < 1, which 
corresponds to to a < 1/2 and £ < 1/2, the correction is 
small. This is the case depicted in Fig. [TJ If £ > 1/2, 
the long- wavelength modes < p\ < ln2\eE\/ir — m 2 
give a large contribution that can only be evaluated in a 
discrete scheme. This case is shown in Fig. [5J The sum 
over the discrete modes provides a delicate cancellation 
between large correlation terms. 

The infinite product ([3]) converges if it is calculated 
over the discrete modes. In order to prove this assertion, 
we proceed with the combinatorial representation of p n : 

1 - 

ax...a n i—1 

The sum here extends over all ordered sets of the n 
pairwise-distinct indices a^. The inclusion in the sum 
of diagonal terms = cx.j increases the value of p n . 
As an upper bound we have the estimate p n < p n = 
po(^L) n /n\, where 

a 

This series converges regardless of whether or not £ 
is small. The cumulative probability appears to be 
bounded: 

00 00 

2_,Pn <^2?n = exp(7L)p - 

n=0 n=0 
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FIG. 3: (color online) The probability distribution (the solid 
curve and empty circles) and the cumulative probability dis- 
tribution (the dashed curve and solid circles) of the number 
of electron-positron pairs, n, produced in a uniform electric 
field of subcritical strength £ = 0.3 for 711 = 10, "/± — 6, and 
7 = 60. Circles 1 show the continuum limit, and curves 2 
correspond to the discrete electron levels. The series of p n 
converge in both cases. The distributions are plotted in the 
log 10 scale. 



FIG. 4: (color online) The probability distribution and the 
cumulative probability distribution of the number of electron- 
positron pairs in a uniform electric field of supercritical 
strength £ = 0.9. The parameters, except for f, and the no- 
tations are the same as in Fig. [3] The series of p n converges 
upon a discrete summation (solid and dashed curves 2). In 
the continuum limit (empty and solid circles 1), the series of 
p n has an asymptotic character. 



The infinite product ([3]) converges; therefore, the p n are 
summed up in the unity, by virtue of Eq. ^ , for any field 
strength. It follows that 1 < exp(7L)po- We remark that 
this inequality is fulfilled in the continuum limit where 
L = L>U 2 (0- 

It is known that limiting procedures modify the an- 
alytic properties of functions. In our case, an infinite 
sequence of zeros in the PGF (Q} has been transmuted 
in a branch cut. The non-analyticity is thus associated 
with the replacement of summation over a by the phase 
space integral. This remark, incidentally, describes how 
to restore analyticity: We have, unsurprisingly, to refuse 
performing calculations in a continuous scheme. This 
prescription is in agreement with the earlier argument 
based on the need to have the cancellation of large cor- 
relation terms, which can be achieved only in a discrete 
scheme. 

One can assume that, in Eq. (HH, n c reveals the max- 
imum value of n for which the calculation of p n in the 
continuum limit is justified for < 1. If this conjecture 
is correct, then expression ([5]) could be valid for arbitrar- 
ily small r^. 

As an illustration, we consider a numerical example of 
the probability distribution both in weak and strong elec- 
tric fields in the continuum limit and for a discrete set of 
levels. The momentum is quantized, as in the nonrela- 
tivistic case for an infinite potential wall with p = im/L, 
where n is a positive integer, L is the size of normalization 
box, and the volume equals V = I? . The boundary con- 
ditions for the Dirac equation have a more complicated 
form. Our goal here is to clarify the principal distinc- 



tion between the continuum limit and a discrete scheme; 
therefore, the relativistic modification of the boundary 
conditions is neglected. 

The probability distributions in the continuum limit 
for £ = 0.3 and £ = 0.9 are shown in Figs. [3] and [4] 
The effective number of degrees of freedom is taken to 
be 7 = 7||7j_, with 711 = \eE\LT/n = 10 and 7j_ = 
\eE\L 2 / {Att 2 ) = 6. In the case £ = 0.3, the distributions 
are smooth, and with the required accuracy the probabil- 
ities sum up to unity. In the second case, where £ = 0.9, 
the series has an asymptotic character. The apparent 
divergence arises starting from n ~ 50. The order-of- 
magnitude estimate (pH) suggests that the correlations 
become large for n > n c ss 15. In the asymptotic region, 
the figure shows the probability for odd n only. The 
probability for even n turns out negative. The violation 
of positive definiteness occurs when the lower-order alter- 
nating terms of the recursion (|19[) become dominant. An 
estimate of the leading alternating term is provided by 
Eq. (|20p . In the region of n greater than 50, a linear semi- 
log dependence is clearly seen; the slope approximately 
equals log(l/r^) ~ 1. We performed the calculation of 
p n in two ways, from the definition (JS|) and via the re- 
cursion (|19p . The results coincide everywhere, including 
the asymptotic region. 

In the asymptotic series, the growing terms are dis- 
carded, starting from the smallest one. The value of the 
first term that is omitted is the intrinsic error of calcu- 
lation. In the situation shown in Fig. |4l the evaluation 
of the series is not numerically satisfactory: Peak values 
of the probability appear inside the region of divergence, 
which makes the intrinsic error large; nor is the cumula- 
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tive probability X)fe<nPfc saturated in the region n < 50. 
With decreasing £, the convergence becomes more consis- 
tent, the asymptotic area moves to the region with high 
n, and, for £ < 1/2, the continuum limit finally provides 
the convergence. 

In discrete schemes, the values of p n are explicitly pos- 
itive. In the continuum limit, this property gets lost. For 
a supercritical field, this fact has been demonstrated in 
Fig. [4] The probability in the tail of the distribution for 
£ < 1/2 is also not positive definite. One can show that 
the values of p n in Fig. [3] are negative, in the continuum 
limit, for even n starting from n = 132. 

The probability distribution is well defined for low n 
in weak fields, where E[n] = 7^ <C n c ~ 7; the variance 
is also small in comparison with n c . In such cases, the 
physically interesting values of n lie in the region below 
n c , where the p n are positive and the continuum limit is 
efficient. 

The results of the discrete summation are also shown 
in Figs. [3] and |U The solid curves correspond to the prob- 
ability distribution, whereas the dashed curves show the 
cumulative probability distribution. The discrete distri- 
butions are smooth everywhere. A strong shift of the 
maximum probability is associated with the replacement 
of the discrete sum by an integral and the relatively low 
value of 7j_ =6. In discrete schemes, the mean value E[n] 
and the variance V[n] can also be expressed in terms of 
the PGF, but because we do not have a closed expression 
for the PGF of a discrete type, this fact does not provide 
any advantages over the direct calculation. Nevertheless, 
the recursion is numerically efficient. The coefficients gi 
of a discrete type are as follows: 

n—l rn—1 v * 

In the continuum limit, we arrive at 

9l = (-l) 1 dxiexpW/t-l)- 1 - 1 , 
Jo 

which coincides with Eq. (fl8|) . 

The long- wavelength part of the spectrum gives a con- 
tribution to gi, which grows exponentially for large I. Ac- 
cordingly, exponential accuracy is required to ensure the 
cancellation of large terms. In the examples considered, 
we carry out calculations with up to 70 significant digits 
in order to cover the range of pair numbers n < 120. In 
the subcritical fields, there is no such difficulty. 

The proximity to the continuum is controlled by the 
value of 7j_. With increasing 7^ and for a constant 7, 
curves 2 in Fig. [4] approach the initial segments of curves 
1, and they remain regular in the region of asymptotic 
divergence of curves 1. 

The probability distribution constructed in Sees. II 
and III corresponds to the one-loop QED effective action 
Q the imaginary part of which determines, in a uniform 
electric field, the pair-production rate. The two-loop 
QED effective action and the corresponding correction 



to pair-production rate have been calculated by Ritus 
[lU , Dittrich and Reuter [l(| and Dittrich and Gies [l!| . 
The effect of the two-loop correction on the probability 
distribution deserves a separate study. 

The generalization of our results to an arbitrary space- 
time dimension d + 1 is straightforward. The PGF takes 
the form of Eq. ([7]), with Li2((l — z)£) being replaced by 
Li s ((l — z)£), where s = (d+l)/2. For double-degenerate 
fermions, the coefficient 7 equals 

2\eE\ s 

ry = VT— - 

7 (2ir) d ' 

where V is the spatial volume. Li s (t) has a branch point 
at t = 1, which leads to the singularity of the PGF at z = 
—r^. The cumulative probability converges for £ < 1/2. 
Therefore, the condition of convergence does not depend 
on the dimension of spacetime. 

Recent interest in Klein's paradox was caused by 
the possibility of an experimental study of the produc- 
tion of quasiparticle pairs using electrostatic barriers 
in graphene [20j. The energy spectrum of quasipar- 
ticles in single-layer graphene has a conical shape, so 
the quasiparticles are described by the Dirac equation 
in 2 + 1 spacetime with massless fermions. In single- 
layer graphene, the branch-point singularity of the PGF 
touches the origin. This gives rise to a strong divergence 
and a failure of the continuum limit approximation. In 
bi-layer graphene, quasiparticles have a gapless parabolic 
spectrum. The analytical properties of the PGF are de- 
termined by the branch cut that also starts at z = 0. 

The Schwinger pair production in multilayer graphene 
is discussed in Refs. |2l| - l23l |. 

The surface of bare strange stars is surrounded by an 
electron layer with a thickness of a few hundred fermi 
[24] ]. The electrons cause a strong electric field in the 
layer. If the temperature of the surface is not zero, 
electron-positron pairs are spontaneously created [25| . 
The effect leads to the radiation of electron-positron 
plasma, which provides a signature for the identification 
of strange stars. On the surfaces, the electric field is 
supercritical, so a discrete scheme could be required for 
quantitative characterization of the pair production in 
strange stars. 

The Schwinger mechanism is used in string models for 
the description of decays of orbitally excited hadrons [26| 
and particle production in high-energy strong interac- 
tions [27]]. The multiplicity of final-state hadrons is de- 
termined by the probability distribution of the number of 
quark-antiquark pairs produced in the color electric field 
of the string. 

V. CONCLUSION 

In this paper, we found an analytical expression for 
the PGF in the continuum limit, and we used it to cal- 
culate the probability distribution, as well as the average 
and variance of the number of electron-positron pairs. 
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The probability for large numbers of pairs was calcu- 
lated using a recursive formula, which we have also pre- 
sented in an explicit analytical form. If the condition 
£ < 1/2 holds, the branch point of the PGF is outside 
the unit circle. The coefficients p n of equation (0, there- 
fore, sum up to unity, confirming the meaning of the PGF 
as a generating function of the probability distribution. 
When £ becomes greater than 1/2, the branch point of 
the PGF steps inside the unit circle, which causes the 
sum of the p n to be divergent. This divergence occurs be- 
cause of an increased contribution of the long- wavelength 
modes, which indicates a failure of the continuum limit 
approximation. The probability in the tail of the distri- 
bution turns out to be not positive definite, for any field 
strength. The unphysical region starts with n ~ n c . In 
the subcritical fields, the cumulative probability is sat- 
urated below the n c , thus justifying results of the con- 
tinuum limit approximation. In the supercritical fields, 
the cumulative probability is not necessarily saturated 
below n c . Whenever this happens, the continuum limit 
approximation fails. 

The PGF, the probability distribution, and the recur- 
sion also were calculated in a discrete scheme. For any 
field strength, the discrete type PGF is an analytic func- 
tion inside the unit circle \z\ < 1; the coefficients p n 



of equation ([5]) are positive for any n; these coefficients 
sum up to unity. In the supercritical fields, discrete sum- 
mation provides, however, only the basic opportunity 
to work with the convergent series and positive-definite 
probability. For quantitative estimates, it is necessary to 
perform calculations with exponential accuracy to guar- 
antee cancellation of large correlation terms. In the sub- 
critical fields, calculations do not require such high accu- 
racy. 

A complete description of the stochastic process is pro- 
vided by a probability distribution. Experimental study 
of the probability distribution of the number of electron- 
positron pairs and/or recursion relations could represent 
the comprehensive test of the Schwinger mechanism. 
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